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ABSTRACT 

The  rate  of  solution  cd  benzoic  acid  in  water  has  been  determined  for  a  system  consist- 
ing of  a  bed  of  acid  spheres  with  water  flowing  through  the  interstices  of  the  bed.  It  is  pro- 
posed to  use  the  results  to  estimate  the  heat  transfer  between  a  sphere  bed  and  a  liquid 
coolant.  The  equation  proposed  for  this  purpose  is 

St  =  2.36  X  10-1  Re-0-3  sc-°-^ 
where  St,  Re  and  Sc  are  the  Stanton,  Reynolds  and  Schmidt  numbers. 
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HEAT  TRANSFER  IN  SPHERE  BEDS 
By  Robert  Bentley,  Frederick  Johnson,  and  Robert  Maurer 

INTRODUCTION 

This  report  represents  an  attempt  to  estimate  the  heat  transfer  properties  of  a  sphere 
bed.  More  specifically  the  problem  is  as  follows.  A  bed  of  metal  or  oxide  spheres,  within 
which  heat  is  being  generated,  is  to  be  cooled  by  water  moving  through  the  interstices  of 
tiie  bed.  For  design  purposes  it  is  necessary  to  know  the  flux  of  heat  from  the  spheres  to 
the  water  as  a  function  of  the  temperature  difference  between  spheres  and  water  and  as 
a  function  of  the  water  velocity  and  bed  geometry    (See  MUC-GY-31).  The  term,  bed  geom- 
etry, includes  such  factors  as  sphere  diameter,  bed  porosity,  and  overall  bed  diameter. 

The  direct  measurement  of  heat  transfer  coefficients  for  such  a  system  is  a  difficult 
and  tedious  task.  Mr.  Young  suggested  that  preliminary  estimates  might  be  made  by  the 
measurement  of  material  transfer  coefficients  and  the  use  of  this  data  to  calculate  heat 
transfer. 

THEORY 

The  heat  transfer  coefficient  h,,is  defined  (see  MUC-GY-28)  by  the  equation 

F„  =hn4T  (1) 

where  Ftj  is  the  flux  of  heat  across  the  solid-liquid  interface,  expressed  in  calories  per 
square  centimeter-second.  AT  is  the  difference  between  the  temperature  of  the  interface 
and  the  average  temperature  of  the  liquid.  A  transfer  velocity  Htj  may  also  be  defined  by 
the  equation 

Fn  =  H„2lC  (2) 

where  4  C  is  the  difference  in  heat  concentration  corresponding  to  the  temperature  dif- 
ference AT.  The  heat  concentration  C  is  defined  by  the  equation 

C  =  cPT  (3) 

and  the  relation  between  hi)   andHu  is 

Hn  =  hfl/Pc  (4) 

The  speciH.c  heat  of  the  liquid  is  c,  and  P  is  its  density.  For  the  analogous  situation  of  a 
material  dissolving  in  a  liquid  stream  from  a  solid  surface,  the  equation 
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Fm  =  Hm4C  (5) 

may  be  written.  Here  Fm  is  the  flux  of  material  across  the  solid- liquid  interface,  expressed 

in  grams  per  square  centimeter-second.  The  difference  between  the  concentration  of  mater- 
ial in  the  liquid  at  the  interface  and  the  average  concentration  of  material  in  the  liquid 
stream  is  ^C.  It  is  assumed  that  the  film  of  liquid  in  contact  with  the  interface  is  saturated 
with  the  material  in  question. 

The  transfer  velocities  for  heat  and  material  defined  by  equations  (2)  and  (5)  wUl  be 
functions  of  the  properties  of  the  system.  If  it  is  assumed  that  the  pertinent  variables  are 
the  viscosity,  density,  and  velocity  of  flow  of  the  liquid,  the  diameter  of  the  spheres,  and 
the  diffusion  coefficient  for  heat  (or  material)  in  the  liquid,  then  a  dimensional  analysis 
indicates  that  the  relationship  between  these  variables  and  the  transfer  velocity  may  be 
written  in  the  following  form. 

St  =  f  (Re,  Sc)  (6) 

The  Stanton  number,  St,  is  (H/V).  Re  is  the  modified  Reynolds  number  (  V  P  d/fi).  Sc  is 
the  Schmidt  number  (/i/PD).  In  the  above  definitions 

V  =  velocity  of  flow  of  the  liquid 

d  =  diameter  of  the  spheres 

H  =  viscosity  of  the  liquid 

D  =  diffusion  coefficient  of  heat  (or  material)  in  the 
liquid 

If  one  works  with  a  given  fluid  and  dissolved  material  at  a  fixed  temperature  the  Schmidt 
number  is  a  constant.  The  diffusion  coefficient  for  heat  in  a  liquid  may  be  related  to  the 
thermal  conductivity  of  the  liquid  by  the  equation. 

D  =  k/Pc  (7) 

The  experimental  problem  is  to  determine  explicitly  equation  (6).  The  assumption  is 
then  made  that  the  relationship  between  the  Stanton,  Reynolds,  and  Schmidt  numbers  which 
is  obtained  experimentally  for  material  transfer  may  be  used  to  calculate  the  heat  transfer 
velocity  if  the  diffusion  coefficient  for  material  is  replaced  by  the  diffusion  coefficient  for 
heat  in  the  Schmidt  number. 

EXPERIMENTAL 

Benzoic  acid  and  w^ter  were  chosen  as  the  materials  to  be  used  in  the  experiment. 
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Benzoic  acid  is  slightly  soluble  in  water  and  the  concentration  is  easily  determined  by 

titration.  The  solubility^  of  benzoic  acid  in  water  has  been  determined  as  a  fimction  of 

temperature  near  room  temperature.  The  dtEfusion  coefficient"  of  benzoic  acid  in  water 
has  been  determined  at ."'  °C. 

Spheres  of  benzoic  acid  were  cast  from  C,P.  material  in  a  brass  mold.  The  best  results 
were  obtained  by  having  the  acid  at  140°C  and  the  mold  at  room  temperature.  The  cast 
spheres  had  a  smooth,  white  surface.  The  average  diameter  of  the  spheres  was  0.635  inches. 

* 

The  form  of  the  sphere  bed  is  shown  in  Figure  1.  The  spheres  were  held  in  a  cylinder 
of  galvanized  iron  closed  at  top  and  bottom  by  wire  screen.  It  was  found  that  the  porosity 
of  the  bed,  which  is  the  fraction  of  the  bed  volume  which  is  void  space,  could  not  be  varied 
but  was  fixed  by  the  diameter  of  the  bed.  Two  4.5  inch  diameter  beds  14  and  21  inches  in 
length  were  examined.  One  8  inch  diameter  bed,  18  inches  in  length,  was  used.  Data  con- 
cerning the  beds  are  given  in  Table  1. 

The  velocity  of  flow  of  the  water  through  the  bed  was  regulated  by  the  head  of  water  and 
the  area  of  the  orifice,  0,  at  the  bottom  of  the  bed.  The  water  velocity,  V,  is  defined  as  the 
quotient  of  the  volume  of  water  passing  through  the  bed  per  unit  time  and  the  cross-sectional 
area  of  the  cylinder  containing  the  bed.  The  volume  of  water  passing  through  the  bed  per 
unit  time  was  measured  by  weighing  the  water  issuing  from  the  bed  in  a  known  interval.  Tap 
water  was  used  in  all  the  experiments. 

Since  the  water  velocity,  V,  does  not  enter  the  Schmidt  number,  a  determination  of  the 
depeiidence  of  the  transfer  velocity,  Hm,  upon  V,  determines  the  relationship  between  the 

Stanton  number  and  the  Reynolds  number.  The  determination  of  Hm  as  a  function  of  V  pro- 
ceeded in  the  following  fashion.  Water,  passed  through  the  bed  for  a  timed  interval,  was 
collected  and  weighed.  From  this  data  the  water  velocity,  V,  was  calculated.  The  concen- 
tration of  benzoic  acid  in  the  water  after  passage  through  the  bed  was  determined  by  titra- 
tion with  0.01  molar  sodium  hydroxide,  using  phenolpthalein  as  an  indicator.  Blank  titrations 
were  made  on  the  tap  water  used  in  the  experiments.  The  total  weight  of  acid  which  dissolved 
into  the  water  stream  was  calculated  from  the  concentration  of  acid  and  the  total  weight  of 
water.  The  flux,  F,  may  then  be  calculated  from  the  above  data,  the  known  total  area  of 
sphere  surface,  and  the  time  of  contact  of  the  water  with  the  bed.  The  concentration  of  dis- 
solved acid  was  in  all  cases  negligible  by  comparison  with  the  solubility  of  the  acid  so  that 
the  concentration  difference  A  C  was  taken  as  the  solubility,  at  the  temperature  of  the  water, 
of  the  acid  in  water. 

A  sample  of  the  data  from  a  typical  e:q>eriment  is  given  in  Table  2.  In  Table  3  is  col- 
lected a  summary  of  all  the  e^qierimental  data.  A  log-log  plot  of  the  transfer  velocity  versus 
the  water  velocity  is  given  in  Figure  2.  From  this  data  one  obtains  the  equation 
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Hnf=7.1  X  10"'*  V^-'^O  (8) 


for  Hrn  and  V. 


In  experiments  of  this  tjrpe  it  is  difficult  to  determine  the  fimctional  dependence  o5  the 
Stanton  number  upon  the  Schmidt  number.  This  is  so  because  none  of  the  quantities  which 
enter  the  Schmidt  number  may  be  readily  varied.  Heat  transfer  experiments  upon  liquids 
flowing  in  pipes  indicate  that  the  Stanton  number  is  inversely  proportional  to  the  Schmidt 

number  raised  to  the  six-tenths  power^.  Experiments  upon  heat  and  material  transfer 
from  solids  to  agitated  liquids  indicate  an  inverse  square  root  dependence  of  the  Stanton 
number  upon  the  Schmidt  number^.  For  flow  of  liquids  normal  to  cylinders  the  Stanton  num- 
ber is  inversely  proportional  to  the  seven-tenths  power  of  the  Schmidt  number-^.  We  have 
assumed  that  the  Stanton  number  is  inversely  proportional  to  the  Schmidt  number  raised  to 
the  six-tenths  power  and  have  rewritten  equation  (8)  in  the  form 

St  =  2.36  x  10"^  Re-0-3  Sc'^-^  (9) 

Since  the  data  for  the  three  beds,  of  varying  height,  diameter,  and  porosity  are  fairly 
well  represented  by  the  single  equation  (8),  there  is  reason  for  belief  that  these  variables 
are  not  of  primary  importance. 

DISCUSSION 

The  assumptions  involved  in  the  use  of  equation  (9)  for  the  calculation  of  heat  transfer 
velocities  must  be  kept  in  mind.  The  first  assimiption  concerns  the  manner  in  which  the 
Schmidt  number  enters  the  equation.  There  is  little  doubt  that  the  inverse  stx-tenttis  power 
is  approximately  correct.  A  better  representation  of  the  behavior  of  the  system  might  be 
obtained,  however,  by  a  change  in  the  exponent  of  the  Schmidt  number  of  ±20%.  Such  a 
change  would  be  of  some  importance  since  the  diffusion  coefficients  for  heat  and  benzoic 
acid  in  water  differ  by  a  factor  of  approximately  one  himdred  (See  MUC-GY-28,  p.3).  In 
passing  by  the  use  of  equation  (9)  from  material  to  heat  transfer,  a  20%  change  in  the  ex- 
ponent of  the  Schmidt  number  changes  the  calculated  heat  transfer  velocity  by  a  factor  of 
approximately  two. 

If  the  results  of  the  present  work  are  compared  with  the  results  of  other  workers,  it 

3 
is  found  that  equation  (9)  predicts  smaller  values  for  the  heat  transfer  velocities.  Furnas 
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and  Saunders  and  Ford^   have  arrived  at  the  following  equations  for  heat  transfer  from 

gases  to  beds  of  broken  solids 

Furnas  St  =  0.8  Re-^-^S  Sc^^-^  (10) 

Saunders  and  Ford  St  =  0.14  Sc"°*^  (11) 

la  each  case,  the  exponent  of  the  Schmidt  number  has  been  assumed.  In  this  connection,  it 
is  worth  noting  that  the  Schmidt  number  for  heat  in  gases  is  of  the  order  of  unity  and  for 
benzoic  acid  in  water  is  of  the  order  of  1000.  This  large  difference  in  the  magnitude  of  the 
Schmidt  nimiber  may  make  a  comparison  of  the  above  equations  with  our  results  impossible. 
The  Schmidt  number  for  heat  in  water  is  approximately  seven.  In  addition  the  data  on  gases 

have  been  taken  for  Reynolds  numbers  of  the  order  of  200.    The  work  of  Gamson^  and  his 
co-workers  indicates  that  this  is  in  the  transition  region  between  streamline  and  turbulent 
flow.  Extrapolation  of  equations  (10)  and  (11)  into  the  region  of  turbulent  flow  is  therefore 
not  permissable.  In  the  present  wort,  the  range  of  Re3molds  number  which  was  covered 
extended  from  200  to  3000. 

Gamson,  working  with  heat  transfer  from  a  sphere  bed  to  a  gas,  found  the  equation 

St  =  1.064  Re"°:'*^  Sc-^-^"^  (12) 


for  Rejmolds  numbers  between  350  and  4000.  The  ejcponent  of  the  Schmidt  number  was 
asstmied.  Equation  (12)  leads  also  to  larger  beat  transfer  velocities  than  our  results  do. 
The  work  of  Gamson  was  performed  upon  very  shallow  beds  which  may  accoimt  for  the 
discrepancies. 
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TABLE  1 

Sphere  Bed  Data 

Bedl 

Bed  2 

Bed  3 

Diameter  of  Sphere  Bed 

11.4 

11.4 

20.3  cm 

Cross -sectional  Area  of  Bed 

102.6 

102.6 

322.1  cm2 

Number  of  Spheres 

1,002. 

1,416. 

4,330. 

Area  of  Sphere  Surface 

8,010. 

11,580. 

35,400.     cm^ 

Height  of  Sphere  Bed 

^5.6 

53.3 

45.7  cm 

Volume  of  Spheres 

2,125. 

3,105. 

9,500.     cm3 

Volume  of  Sphere  Bed 

3,647. 

5,460. 

14,720.     cm^ 

Porosity 

43.% 

43.  % 

36.  % 
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TABLE  3 

Summary  of  Data 

* 

FxlO^ 

,      Z\C  X  10^ 

H    x.lO^ 

V 

I     ^               \ 

cm'' 

,    cm 

^     sec      ' 

sec 

^        2                 ' 
cm''  sec 

Bedl 

1.69 

2.69 

6.3 

1.07 

2.49 

2.46 

10.1 

1.95 

6.86 

2.47 

27.8 

4.50 

8.53 

2.95 

28.9 

6.75 

8.79 

3.01 

29.2 

7.95 

11.97 

2.64 

45.3 

14.2 

16.44 

2.52 

65.3 

21.1 

Bed  2 

2.11 

2.45 

8.6 

1.25 

5.44 

2.47 

22.0 

4.60 

6.57 

2.50 

26.3 

6.71 

10.6 

2.46 

43.1 

13.6 

12.8 

2.44 

52.5 

15.4 

Beds 

1.87 

2.46 

7.6 

1.07 

3.05 

2.42 

12.6 

2.32 

4.95 

2.46 

20.1 

4.44 

6.93 

2.45 

28.3 

8.44 

8.85 

2.45 

35.1 

14.7 

11.8 

2.59 

45.7 

17:3 

Diffusion  coefficient  of  benzoic  acid  in  water 


cm' 


at  250c      (6) 


D  =  1.1  X  10-5 

sec 

D  =  8.3  X  10-6  cm^    at  IS^C     (Calculated;  see  CP-2883) 
sec 

H  =  1.14  X  10-^2  poige  3t  i^Oq 

Calculation  of  Reynolds  and  Schmidt  numbers  peiiormed  for 
a  temperature  of  15°C. 
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APPENDIX 


Heat  transfer  to  liquids  flowing  in  smooth  walled  pipes  of  circular  cross  section  has 
been  extensively  investigated  and  for  the  region  of  well  developed  turbulent  flow  (Re>  10,000) 
the  following  equation  is  well  established. 

St  =  2.30  X  10-2  Re"°-2  Sc"°-^  (13) 

Material  transfer  velocities  were  measured  in  a  circular  pipe  as  a  function  of  the  linear 
water  velocity  to  determine  if  an  equation  similar  to  (13)  existed.  Figure  3  (a)  is  a  section 
of  the  benzoic  acid  pipe  used  in  this  experiment.  A  brass  core  was  centered  in  a  brass  tube 
by  means  of  end  caps  and  benzoic  acid  cast  in  the  space  between  the  core  and  the  outer 
brass  shell.  The  brass  core  was  then  withdrawn.  Two  such  sections  were  used  in  the  ar- 
rangement illustrated  in  Figure  3  (b).  The  water  was  led  from  a  reservoir  to  the  benzoic 
acid  pipes  by  means  of  a  brass  pipe  approximately  six  feet  long  and  of  the  same  internal 
diameter  as  the  acid  pipe  units.  The  details  of  procedure  for  the  determination  of  the 
transfer  velocity  were  the  same  as  in  the  sphere  bed  experiments.  Table  4  is  a  simnmary 
of  the  data.  Solution  of  the  benzoic  acid  in  the  water  during  the  experiment  increased  ap- 
preciably the  diameter  of  the  pipe  section.  Assuming  that  the  change  in  diameter  was 
constant  for  each  run,  a  corrected  diameter  was  used  in  computing  the  Reynolds  number 
and  surface  area  for  each  e:g)eriment.  In  Table  4  are  given  the  initial  and  final  areas  of 
the  benzoic  acid  surface.  Figure  4  is  a  log-log  plot  of  the  material  transfer  velocity  versus 
the  linear  water  velocity.  The  straight  line  of  Figure  4  represents  the  equation 

H  =  1.23  X  10"'*  V°-^  (14) 

Making  the  usual  assumption  concerning  the  dependence  of  the  Stanton  number  upon  the 
Schmidt  number,  equation  (14)  may  be  rewritten  as 

St  =  2.04  X  10-2  Re-0-2  Sc'^-^  (15) 

Comparison  of  this  result  with  equation  (13)  for  heat  transfer  indicates  satisfactory  agree- 
ment. This  result  lends  considerable  weight  to  the  proposed  use  of  material  transfer  data 
for  the  estimation  of  heat  transfer  performance. 

It  is  rather  surprising  that  the  data  of  Figure  4  fit  a  straight  line  plot  of  slope  0.8  as 
well  as  they  do.  For  Reynolds  numbers  less  than  10,000  the  heat  transfer  data  indicate  a 
greater  slope.  A  partial  explanation  may  lie  in  the  fact  that  in  the  course  of  solution  of 
the  acid  the  surface  was  appreciably  roughened.  The  excellent  agreement  of  equations  (13) 
and  (15)  must  therefore  be  considered  somewhat  fortuitous. 
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TABLE  4 
Summary  of  Results  for  Pipe  Experiment 

Diameter  of  .pipe:  0.95  cm 
Cross  sectional  area  of  pipe:  .709  cm^ 
Total  Length  of  acid  sections:  58.5  cm 
Initial  area  of  acid:  174.2  cm^ 
Final  area  of  acid:  237.1  cm^ 


xlO^ 

h ) 

m^  sec 

Ac  X  103 

(     ^             ) 
cm3 

HxlO^ 

(      ""^         ) 
^     sec 

V 

(     ^^    ) 

H.  -  ( 7') 

^     sec   > 

3.89 

2.50 

15.55 

25.2 

2,390 

6.50 

2.48 

26.20 

35.3 

3,350 

7.70 

2.48 

31.01 

55.2 

5,260 

3.28 

2.49 

53.20 

123.3 

11,600 

8.30 

2.50 

73.10 

188.9 

17,900 

8.10 

2.55 

110.0  • 

270 

25,700 

8.70 

2.53 

113.4 

295 

28,000 

;1.90 

2,53 

86.60 

198 

18,800 

.8.25 

2.53 

72.10 

147.9 

14,100 

.1.00 

2.55 

43.10 

78 

7,410 

7.45 

2.56 

29.10 

51.9 

4,930 

3.98 

2.58 

15.41 

36.9 

3,500 

-  14  ^- 
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